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Abstract 

The microbiota of four individual deep water sponges, Lissodendoryx diversichela, Poecillastra compressa, Inflatella pellicula, 
and Stelletta normani, together with surrounding seawater were analysed by pyrosequencing of a region of the 16S rRNA 
gene common to Bacteria and Archaea. Due to sampling constraints at depths below 700 m duplicate samples were not 
collected. The microbial communities of L diversichela, P. compressa and /. pellicula were typical of low microbial abundance 
(LMA) sponges while S. normani had a community more typical of high microbial abundance (HMA) sponges. Analysis of the 
deep sea sponge microbiota revealed that the three LMA-like sponges shared a set of abundant OTUs that were distinct 
from those associated with sponges from shallow waters. Comparison of the pyrosequencing data with that from shallow 
water sponges revealed that the microbial communities of all sponges analysed have similar archaea! populations but that 
the bacterial populations of the deep sea sponges were distinct. Further analysis of the common and abundant OTUs from 
the three LMA-like sponges placed them within the groups of ammonia oxidising Archaea {Thaumarchaeota) and sulphur 
oxidising y-Proteobacteria (Chromatiales). Reads from these two groups made up over 70% of all 16S rRNA genes detected 
from the three LMA-like sponge samples, providing evidence of a putative common microbial assemblage associated with 
deep sea LMA sponges. 
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Introduction 

Marine sponges are important members of marine benthic 
communities throughout polar, tropical and temperate oceans. 
They are sessile filter feeding animals, capable of filtering up to 50 
thousand litres of seawater per litre of sponge per day [1]. Marine 
sponges have gained much attention during recent years due to 
their remarkably dense and diverse community of bacterial, 
archaeal and eukaryotic microorganisms. This microbiota con- 
tributes to sponge biology in many ways, such as providing a 
chemical defence mechanism, carbon and nitrogen cycling and as 
a food source [2] [3]. Sponges can also be grouped according to 
the density of bacteria within their tissues into high microbial 
abundance (HMA) sponges and low microbial abundance (LMA) 
sponges [4] with HMA sponges reported to have microbial 
densities of lO'^-lO"' bacteria per g of tissue while LMA sponges 
have 10'- 10 bacteria per g of tissue. Great efforts have been 
made to characterise the diversity of the microbial assemblages in 
shallow water sponges, and in HMA sponges in particular, using 
both culture dependent [5-8] and culture independent approaches 
[9-11]. Up to 2007, 15 bacterial phyla (including the candidate 
phylum Poribacteria), 2 major archaeal lineages and many microbial 
eukaryotes had been reported from marine sponges [2]. In recent 
years the application of "next generation" sequencing has allowed 



access to the so called "rare biosphere" [12] and increased the 
number of bacterial phyla detected in sponges to more than 30. 
Marine sponges from the Great Barrier Reef [3], the Red Sea 
[13], the Mediterranean [14], the northern Atiantic [15,16], the 
Caribbean [17], Brazil [18] and worldwide [19] have been studied 
for their microbial diversity. More interesting than the sheer 
diversity of microbial communities which have been found in 
sponges are the presence of sponge-specific microorganisms, i.e. 
OTUs found almost exclusively in sponges [20] . Li a comprehen- 
sive study of sponge-microbe associations, Schmitt et al. have 
distinguished between core, variable and species-specific assem- 
blages in sponges. Interestingly, only a very small proportion of 
90%, 95% and 97% OTUs was shared between different sponge 
species [19]. This "core" community of microbial OTUs found in 
most studied sponges implies a horizontal transfer of sponge- 
associated microbial diversity through the surrounding seawater. 
In previous studies, by comparison of larvae and adult sponges, 
evidence for vertical symbiont transfer has also been shown [3]. 
Thus potentially both vertical and horizontal transfer is involved in 
shaping sponge-associated microbial communities. Other studies 
have focused on describing the community structures in diverse 
sponges [3,15], including archaeal diversity [13], seasonal 
variations in the community structure [17] and functional analysis 
of the sponge metagenome [16,18]. 
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Relative to their shallow water counterparts, little is known 
about the microbial diversity of deep water sponges, due in no 
small part to the inherent technical difficulties in obtaining 
specimens from below 1 000 m. Nevertheless, in culture dependent 
studies Romanenko and co-workers reported the isolation of two 
new bacterial species from a deep sea sponge [21,22], while Brtick 
and colleagues identified an Entotheonella species in Discodermia 
dissoluta from a depth of 150 m [23] and later characterized the 
culturable anaerobes from Geodia sp. samples from depths of 
~200-350 m [24]. A culture independent study on Polymastia cf. 
corticata sampled at a depth of ~ 1 1 00 m revealed that bacteria 
previously found in shallow water sponges are also present in deep 
water sponges and that the bacterial community has a spatial 
distribution in the sponge [25]; a phenomenon which has also 
been described for a shallow water sponge [26]. Deep water 
sponges, including a Lissodendoryx sp., [27,28] have also proven 
fruitful sources of novel bioactive compounds, many of which are 
likely to be of microbial origin. 

The potential roles of most microorganisms within the sponge 
microbiota are, as yet, largely unknown. However recent studies 
using biochemical and metatranscriptomic approaches strongly 
suggest that ammonia-oxidising Archaea (AOA) are actively 
involved in nitrification within sponges [16,29]. Roles for other 
members of the sponge microbiota in processes such as sulphur 
oxidation and the provision of chemical defence systems are also 
likely but currently there is little direct evidence of these. 

In order to increase our understanding of the microbial 
communities associated with deep water marine sponges and to 
assess potential similarities between the microbiota of deep water 
species, this study has applied the pyrosequencing approach to 
analyse the microbiome of sponge samples from the bathypelagic 
zone. Four northern deep water species, /. pellicula, S. normani, L. 
diversichela and P. compressa which have not previously been studied 
for their associated microbiota, were collected from deep water 
canyon systems in the North Atiantic at deptiis of 748, 1 350, 1 350 
and 1469 m respectively. The primer pair targeting the V5-V6 
region of 16S rRNA genes which is common between both Archaea 
and Bacteria was chosen because Archaea have been reported to be 
particularly abundant in deep water marine sponges [30]. Also, 
this primer pair yields sequence lengths of about 280 bp which 
enables the classification of sequence reads into lower taxonomic 
levels. The results of these analyses have led to the characterisation 
of the microbial communities associated with these sponge samples 
and the comparison of these communities has revealed a putative 
common microbial assemblage associated with deep sea LMA-like 
sponges. A fuller understanding of the deep sea sponge microbial 
community is a first step to understanding the roles these microbes 
play in this as yet poorly understood environment. 

Results 

Sponge sampling 

Four sponge samples and water samples were obtained from 
depths of 700 m to 1500 m (Table 1). Sponge samples (Figure 1) 
were identified as L. diversichela (LD), S. normani (SN), P. compressa 
(PC) and /. pellicula (IP). Duplicate samples were not available for 
analyses due to the difficulties associated with sampling at such 
depths. The cortex and the choanosome of S. normani were 
processed as two separate samples in order to characterize any 
potential spatial distribution of bacteria in this sponge although for 
most analyses the data from these samples were combined. The 
outer layer of the sponge was cleaned carefully with a sterile 
scalpel in order to remove any sediment attached to the sponge. 



Deep Sea Marine Sponge Microbiome 




Figure 1. Sponge samples L. diversichela (LD) and 5. normani 
(SN). 

doi:1 0.1 371/journal.pone.0091 092.g001 

Taxonomic richness 

A total number of 81,244 individual 16S rRNA sequence reads 
were obtained by pyrosequencing. Errors in the dataset were 
analysed and corrected using Acacia [31]. The average sequence 
length of all quality filtered sequences was ~280 bp. Following 
these quality checks a total of 70,582 sequences were analysed 
(Table 1) using QIIME [32]. Operational taxonomic units (OTUs) 
were selected at 97%, 94%, and 90% similarities and potentially 
chimeric sequences were removed using ChimeraSlayer following 
aUgnment. A total of 6,357 (97% similarity) OTUs were obtained 
from the sponge and water samples with 4,504 OTUs from the 
sponge samples alone (Table 2). The sample from sponge Stelleta 
normani had the most diverse population with 3,942 97% OTUs 
detected in this single sponge sample. The other sponge samples 
had OTUs ranging from 172 to 247, while OTUs in the water 
samples ranged from 382 to 1,213. 

The rank abundance curves at 97% sequence similarity (Figure 
SI) show that the microbial communities of the samples from 
sponges L. diversichela and /. pellicula are dominated by a relatively 
small number of OTUs, demonstrated through the steep slope. 
The communities of the other samples are more evenly distributed 
with S. normani having the most even distribution. The rarefaction 
curve (Figure S2) of the microbial associates of the S. normani 
sample does not reach a plateau, indicating that even with over 
22,000 reads there is an undersampling of the microbial 
biodiversity. The rarefaction curves of the other samples indicated 
that a greater proportion of the biodiversity was sampled although 
two of the water samples also indicated that there was some under 
sampling of these. 

Bacterial vs. archaeal diversity 

The relative abundance of archaeal reads in the sponge samples 
ranged from 4% in LD, 19% in SN, 48% in IP and up to 65% in 
PC, while in the accompanying surrounding seawater samples the 
Archaea made up 36-38% of all reads. Archaea were more often 
found in the cortex of the S. normani sample (28% relative 
abundance) than in the choanosome (10%). In all samples, 
bacterial reads made up the remaining proportion with a 
negligible amount of sequencing reads not classified into either 
of the two prokaryotic domains of life. 

Archaeal diversity 

In all sponge and water samples the two archaeal phyla 
Emyarchaeota and Thaumarchaeota were present (Figure 2). Crenarch- 
aeota were detected at very low levels (two single reads) in two of 
the samples (W2 and W3). A total of 192 (97%) OTUs were 
classified as Archaea. Among the sponge samples there were 61 
(97%) OTUs within the Archaea with the Thaumarchaeota being the 
most diverse Archaea present (46 OTUs). Over 95% of all archaeal 
reads in the sponge samples were Thaumarchaeota with the 
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Table 1. Sponges and water samples and pyrosequence data. 





Sample 




latitude 


longitude 


depth [m] 


initial reads 


Filtered reads* 


/, pellicula 


IP 


54.0015 


-12.3100 


748 


10070 


8470 


L diversichela 


LD 


54.0584 


-12.5469 


1350 


6423 


6093 


P. compressa 


PC 


54.0633 


-12.4131 


1469 


4445 


1976 


S. normani 


5N 


54.0613 


-12.5518 


1350 


26668 


22475 


Seawater 


W-1 


54.0015 


-12.3100 


748 


6207 


5863 


Seawater 


W-2 


54.0630 


-12.4165 


1500 


18328 


17044 


Seawater 


W-3 


54.0584 


-12.5469 


1350 


9103 


8661 



The sponge and water samples were collected at depths from 700 m to 1500 m in canyon regions north of Porcupine Bank in the N. Atlantic. 
*After applying Chimera Slayer from reads grouped into OTUs at 97% similarity. 
doi:1 0.1 371 /journal.pone.0091 092.t001 



remainder classified as Euryarchaeota. Tlie seawater samples 
contained both Thaumarchaeota and Euryarchaeota with the Thau- 
marchaeota fi-action consisting of 17% (W-1), 62% (W-2) and 68% 
(W-3) of archaeal reads. Of the Thaumarchaeota fraction the vast 
majority of the seawater-derived reads and those from sponge IP 
were classified in the candidate genus Mtrosopumilus, whereas reads 
from the remaining sponges were mainly classified in the family 
Cenarchaeaceae. In each sponge sample over 70% of the Thaumarch- 
aeota fraction was made up of up to three (97%) OTUs. Among 
reads classified as Euryarchaeota there were no abundant OTUs 
from the sponge samples, with the water samples each containing 
several OTUs in the marine group II and III Thermoplasmata. 

To allow further analysis of the archaeal members of the sponge 
associated communities, longer fragments of the archaeal 16S 
rRNA genes were amplified and cloned. Each of these sponge 
hbraries was found to contain clones that were close matches to 
abundant amplicons from the respective pyrosequencing libraries. 
AU the archaeal 16S rRNA gene sequences recovered from the 
clone libraries were analysed by BLAST analysis and by 
taxonomic classification with Greengenes, RDP and Silva 
pipelines. Phylogenetic analysis (Figure 3) and Greengenes 
taxonomy (Table 3) showed that all sponge archaeal clones 
grouped with the Marine Group I Thaumarchaeota. Phylogenetic 
analysis revealed two distinct clades; a clade of sponge-derived 
reads from sponges /. pellicula, P. compressa and S. normani that 
clustered together with Mtrosopumilus maritimus and the sponge 
symbiont Cenarchaeum symbiosum; and a clade consisting of clones 
from sponges L. diversichela and P. compressa that clustered with 
seawater-derived archaeal reads (Figure 3). 

Bacterial diversity 

A total of 6, 160 bacterial OTUs were detected at 97% similarity 
with 4,506 bacterial OTUs (97%) present within the sponge reads 
(Figure 4). These OTUs were assigned to a total of 35 bacterial 
phyla (including candidate phyla) with 28 phyla detected in the 
sponge-derived reads. The most diverse phylum of bacteria within 
the sponge samples were the Chlorojlexi with 1,487 individual 
OTUs (97%) present (this was reduced to 118 OTUs when 
sequences were grouped at 90% sequence similarity), the 
Proteobacteria were the next most diverse phylum with 722 OTUs 
(97%). The proteobacterial OTUs could be further subdivided 
into 317 y-Proteobacteria, 128 ^-Proteobacteria and 167 a-Proteobacteria. 
Other diverse groups of bacteria present included the Acidobacteria 
with 321 OTUs (97%), the Poribacteria with 1 14 OTUs (97%), the 
Actinobacteria with 110 OTUs (97%), Gemmatimonadetes with 54 
OTUs (97%) and Bacteroidetes with 48 OTUs (97%). A large 
number of 97% OTUs (1297) were classified only at the domain 



level but these were conflated to 234 when a 90% similarity cut off 
was used. A summary of OTU numbers at different similarity 
groupings is shown in Table 2. 

A number of candidate divisions /Phyla were also detected in 
addition to the Poribacteria; GN02, OP3, SAR406, TM7, and ZB3 
were all detected with low read abundance and diversity, with 
candidate divisions/Phyla SBR1093 and WS3 being more diverse 
(23 OTUs (97%) and 25 OTUs (97%) respectively) in the sponge 
samples. 

Analysis at lower taxonomic levels 

Analysis at lower taxonomic levels revealed that the Thaumarch- 
aea present in all sponges were classified in the genera Cenarchaeum 
and Mtrosopumilus (Figure 3). From genomic analysis of the sponge 
symbiont Cenarchaeum symbiosum and physiological studies of 
cultured Mtrosopumilus these species are known to belong to the 
group of ammonia-oxidising archaea (AO A) capable of oxidising 
ammonia to nitrite [33-35]. Other microorganisms present within 
the sponges with potential contribution to nitrogen cycling are the 
nitrite oxidising bacteria Mitrospiraceae, Mitrospina and the ammonia 
oxidising bacteria Mtrosomonadales. Reads classified to Mtrospira 
were found in three of the sponge samples though they were more 
prevalent in SN and IP, whereas Mitrospina classified reads were 
found in all sponge samples. 

Of the other common phyla detected within the sponges 
Acidobacteria group 6, and Chlorojlexi SAR202 were found in all 
sponge samples and were especially common in SN (Figure 3). 
These groups are typical members of the sponge microbiota and 
have been found in numerous studies of sponges from diverse 
habitats including sponges from shallow, tropical environments 
[13] and from deeper, cold water environments [16]. However all 
sponge samples contained OTUs consisting of large numbers of 
reads that were classified only at the level of the class y- 
Proteobacteria. When the data was analysed at 90% similarity 
grouping, many of these OTUs grouped together to form single 
OTUs. 

As these unclassified y-proteobacterial reads appeared to 
represent an abundant OTU within the sponges we sought to 
study these in greater detail. To allow for more robust 
phylogenetic analyses than are possible with the shorter reads 
generated by pyrosequencing near full length 16S rRNA libraries 
were prepared from each sponge metagenome with several clones 
(4—30) sequenced from each. Comparison to the pyrosequencing 
data showed that each of these small libraries contained 
representative clones of the most common OTUs detected in the 
pyrosequencing data (including unclassified y-Proteobacteria, Mitros- 
pira, Chlorojlexi and others). These sequences were analysed further 



PLOS ONE I www.plosone.org 



3 



March 2014 | Volume 9 | Issue 3 | e91092 



Deep Sea Marine Sponge Microbiome 

by BLAST analysis, by taxonomic assignment using RDP, 
Greengenes and Silva pipelines and by alignment and phyloge- 
netic analyses. Greengenes taxonomic assignment classified the 
previously unclassified y-proteobacterial OTUs within the Chro- 
matiales or with sulphur oxidising symbionts (Table 3). Likewise 
phylogenetic analysis of these clones showed that all these 16S 
rRNA sequences were found to cluster within the Chromatiales 
group of y-Proteobacteria together with other sponge-derived 
sequences and endosymbionts from tubeworms (Figure 5). AU 
cultured members of these groups and characterised symbionts are 
believed to be involved in sulphide oxidation [36-38]. Other 
sequenced clones confirmed the presence of multiple Chlorojlexi 
from S. normani, Mtrospira and JVitrospina bacteria together with 
other groups detected from the pyrosequencing analysis [Acidimi- 
crobiales, Gemmatimonadetes and Acidobackria). 

Analyses of sponge communities 

Robust comparison of sponge species specific microbial 
communities was not possible due to the lack of duplicate sponge 
samples, however shared features of the deep sea sponge microbial 
communities were examined. In order to determine whether there 
was a shared deep-sea sponge microbial community the compo- 
sitions of the sponge microbial communities were analysed by 
direct comparison of shared OTUs. Classified reads were first 
grouped into OTUs at 97%, 94% and 90% sequence similarity, 
OTUs shared between sponge samples were then analysed to 
determine if there was a common microbiota. For the purposes of 
this study the shared microbiota was defined as OTUs that were 
present in at least 3 of the 4 sponge species analysed (Table 4). 

At all sequence similarity groupings (90%, 94% and 97%) over 
70% of all OTUs were present within just a single sponge sample. 
At 97% similarity grouping there were very few shared OTUs with 
only 3 OTUs shared between 3 or more sponge samples, and 25 
present in 2 or more. The analysis of shared OTUs at 94% 
similarity gave 23 OTUs that were present in 3 or more sponge 
samples and at 90% similarity there were 37 OTUs present in 3 or 
more sponge samples (Table 4). While the number of OTUs 
shared between the sponge samples was low compared to the 
overall number of OTUs present, when the number of reads 
associated with the shared OTUs were analysed it was found that 
the common microbiota (at 90% OTU similarity) comprised 25% 
(SN), 86% (IP), 87% (PC), and 95% (LD) of the total number of 
analysed reads for each sponge sample. The microbiota of the 
sponge and seawater samples was also compared using UniFrac 
UPGMA cluster analysis. The relationship between the samples 
shows that the three water samples from the different depths and 
locations have similar microbial communities and form a cluster 
(Figure 6). The sponge samples clustered into two groups, one 
consisting of the choanosome and cortex from the two samples of 
SN and the other cluster consisting of the microbiota of IP, PC 
and LD. Thus the microbial communities of the two tissue types 
from the S. normani sample are highly similar as are the 
communities from the single samples of sponges /. pellicula, P. 
compressa, and L. diversichela (Figure 6). 

The shared microbiota was found to contain 3 OTUs (90%) 
classified as Thaumarchaeota (also in Figure 3) comprising a total of 
over 95% of all reads characterised as Archaea. In addition the 
shared microbiota was found to contain an OTU within the group 
oiy~Proteobacteria described above as grouping with the Chromatiales. 
This OTU comprised 69%, 41% and 90% of the bacterial reads 
from IP, PC and LD respectively. A related OTU within the 
Chromatiales group (identified by BLAST and phylogenetic analysis) 
from SN comprised 6% of the bacterial reads from this sponge 
sample (Figure 5). While the common sponge community 
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Figure 2. Archaeal diversity in sponge and seawater samples. OTUs were grouped at 97% similarity and taxonomy assigned using the RDP 
classifier with the Greengenes database of assigned sequences. Refer to Table 1 for sample abbreviation. 
doi:1 0.1 371/journal.pone.0091 092.g002 
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uncultured Nitrospira sp. Xestospongia testudinaria (HQ270207) 
IPE7 

Nitrospira marina str. 295 (X82559) 

Chloroflexi (including 4 SN clones) 
Archaea 



Chromatiales 



y-Proteobacteria 



Nitrospira 



Figure 3. PKiyiogenetic relationship of archaeal 16S rRNA clones from deep sea sponges. 16S rRNA sequences were determined from 
cloned PCR amplicons (also see Table 3). Sequences were aligned using PyNast and neighbour joining phylogenetic tree was constructed using 
MEGA5. Bootstrap values (500 replicates) greater than 50% are shown next to the branches. 
doi:1 0.1 371/journal.pone.0091 092.g003 
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■ SAR406; AB16 

□ SBR1093; EC214 

■ gammaproteobacteria; other 
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Figure 4. Bacterial diversity in sponge and seawater samples. OTUs were grouped at 97% similarity and taxonomy assigned using the RDP 
classifier with the Greengenes database of assigned sequences. Refer to Table 1 for sample abbreviation. 
doi:1 0.1 371 /journal.pone.0091 092.g004 



described above accounted for over 75% of all reads in sponge 
samples from IP, PC, and LD, the SN sponge sample had a more 
complex community with the presence of many reads classified as 
Acidobacteria, Poribacteria, Chioroflexi, Mitrospira and candidate 
SBR1093 (Figure 3). The high diversity and overall makeup of 
the microbial community of the SN sample is similar to that found 
in high microbial abundance (HMA) sponges from other marine 
environments [39] and this data clearly supports the classification 
of SN as a deep water HMA sponge. In contrast, the dominance of 
IP, PC and LD by a few OTUs is typical of low microbial 
abundance (LMA) sponges [40] and these three sponge samples 
are likely to be examples of deep sea LMA sponges. 

Comparison to other sponge communities 

To compare the microbial communities of the deep sea sponges 
that were the subject of this study with other sponge microbial 
communities the dataset from a study of sponges collected from 
the Red Sea at depths of 8-19 m was analysed [13]. These data 
were generated using the same primer pairs as used in this study so 
a direct comparison of the data was possible. The Red Sea sponge 
data were reanalysed together with the deep sea data. 

Comparison of shared OTUs grouped at 90% similarity 
indicated that there were similarities in the archaeal communities 
(Table 5). One archaeal OTU, classified as C. sjmhiosum was 
present in all samples making up a large proportion (20-98%) of 
all reads classified as Archaea, while two other OTUs classified as 
Archaea were detected in 6 of the 7 sponge samples. These three 
common archaeal OTUs together accounted for 65-99% of total 
archaeal reads (Table 5A). One bacterial OTU, classified as 
Acidobacteria- 6, was found in all samples, however the abundance of 
these reads was low (<0. 1 %) in three of the sponge samples. Five 
other bacterial OTUs were found in 6 of the 7 sponge samples. 



however all these shared bacterial OTUs made up only 1—7% of 
the total reads (Table 5A). In contrast, analysis of the Red Sea 
shallow water sponges revealed 103 shared OTUs in total, 6 
archaeal OTUs and 97 bacterial OTUs, together these shared 
OTUs constituted 82-93% of the archaeal reads and 63-79% of 
the bacterial reads from the Red Sea shallow water sponges 
indicating a high degree of similarity in the microbial communities 
of these sponges (Table 5B). Of the 97 bacterial OTUs common to 
the Red Sea shallow water sponges, only 7- 1 1 were present in the 
deep sea LMA sponges and these made up 1-23% of the total 
bacterial reads (Table 5B). A total of 27 OTUs were found to be 
common to the deep sea LMA sponges, 5 archaeal and 22 
bacterial. These shared bacterial OTUs made up a total of 67— 
94% of the deep sea bacterial reads (Table 5C). Of the bacterial 
OTUs common to the deep sea LMA sponges, 3-14 were also 
detected in the Red Sea shallow water sponges, however these 
were found to make up only a small component (2-5%) of the 
bacterial population of the Red Sea sponges, highlighting 
differences in the microbial community structures. 

The microbial communities of deep sea and Red Sea shallow 
water sponges were also analysed using UniFrac UPGMA cluster 
analysis (Figure 7). The results of this largely supported the analysis 
of the shared OTUs, the microbiota of the three LMA-type deep 
water sponges (LD, IP and PC) were found to be more similar to 
each other than to other sponge or water samples with the 
microbiota of the LMA Red Sea sponge Sfylissa carteri being more 
closely related than other sponge or water samples. The microbiota 
of SN was found to be more distinct but again was more related to 
the microbiota of the Red Sea sponges Xestospongia testudinaria (a 
confirmed HMA sponge) and Hjrtios erectus (a presumed HMA 
sponge). The UniFrac analysis thus supports the classification of the 
sponges LD, IP and PC as LMA and SN as HMA. 
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Figure 5. Phylogenetic relationship of bacterial 16S rRNA clones from deep sea sponges. 16S rRNA sequences were determined from 
cloned PCR amplicons (also see Table 3). Sequences were aligned using PyNast and neighbour joining phylogenetic tree was constructed using 
MEGA5. One representative sequence per OTU was used for alignment and construction of tree. Bootstrap values (500 replicates) greater than 50% 
are shown next to the branches. 
doi:1 0.1 371/journal.pone.0091 092.g005 



Discussion 

Marine sponges have long been known to harbour a wide range 
of microbes and while the microbial community of shallow water 
sponge.s has been studied extensively, relatively little is known 
about the microbes associated with deep water sponges, especially 
from the bathypelagic and deeper zones. By applying the 454 
pyrosequencing approach to assess the microbial communities of 
deep-water sponge samples this study sheds further light on 
microbes associated with these sponges. 

Bacterial community 

The lack of duplicate sponge samples in this study means that 
robust conclusions about the microbial communities associated 
with individual deep seawater sponge species are not feasible. 
Howewer the shared microbiota of these deep sea sponge species 
has been analysed to determine whether there is a specilic 
microbial community associated with deep water sponges. The 
analysis of shared OTUs between the deep sea sponge samples 
revealed that most of the OTUs present were unique to individual 
sponge samples, however there was also a number of OTUs that 
were shared between at least 3 of the 4 deep water sponge species 
investigated, suggesting a shared core community in these deep 



water sponges. Although sequence abundance can not be 
considered truly quantitative due to biases in PCR amplification 
and differences in 16S rRNA gene copy number these OTUs were 
found, in some cases, to constitute very large proportions of the 
total reads in the study, indicating that they are likely to make up 
substantial parts of the sponge microbial communities. Analysis of 
pyrosequencing read abundance has proved insightful in a recent 
study describing and demonstrating differences in microbial-host 
specificity in LMA sponges [41]. Microbes within these abundant 
groups included the AOA Thaumarchaeota, a group which contains 
the sponge symbiont C. symbiosum. A second group of microbes 
present in all sponge samples were identified as y-Proteobackria 
within the Chromatiale.s group. In three of the sponge samples IP, 
LD and PC these two microbial groupings made up over 70% of 
all reads. This lower diversity microbial community is typical of 
low-microbial abundance (LMA) sponges, however other studies of 
LMA sponges have concluded that each sponge species harbours a 
unique bacterial community [40], while the data here suggest a 
potential common microbiota for deep sea LMA sponges. The 
data from this study also confirm other recent reports that LMA 
sponges have microbial communities that are distinct from and are 
not merely a reflection of surrounding seawater [40,41]. Direct 
comparison of pyrosequencing data, including shallow water LMA 



PLOS ONE I www.plosone.org 



8 



March 2014 | Volume 9 | Issue 3 | e91092 



Deep Sea Marine Sponge Microbiome 



SNCo 

SNCh 

-PC 
— IP 



LD 

- W-2 

- W-3 



- W-1 



3 ^ 
° 

O -D 



Figure 6. Relationship between microbial communities of deep 
sea sponges and seawater by UniFrac UPGMA clustering. See 

table 1 for sample abbreviations. The pyrosequencing data for the 
Stelletta normani cortex (SNCo) and choanosome (SNCh) samples were 
analysed separately. Jackknife values are shown at nodes. The scale bar 
indicates the Unifrac distance. 
doi:10.1371/journal.pone.0091092.g006 

sponges, with this dataset revealed that the microbial communities 
of all the sponges have similar archaeal populations but that the 
bacterial populations of the deep sea sponges are distinct. Other 
non-pyrosequencing based studies of LMA sponges collected from 
shallow temperate and tropical waters show that LMA sponges 
typically have a single large bacterial OTU that can be 
cyanobacterial or proteobacterial, and that appears to be species 
specific [40]. This is in contrast to the deep sea LMA sponges 
which, from the data presented here, appear to have a more 
conserved bacterial population structure. Two additional studies of 
deep sea sponge microbiota provide additional support for this; 
analysis of the microbial community of a single sample of a deep 
sea Polymastia sponge by DGGE analysis of 16S rRNA genes and 
aprA (adenylyl-sulphate reductase) genes revealed the presence of 
sulphur-oxidising y-Proteobacteria (Chromatiaks) and marine group I 
Thaumarchaeota [25]. Likewise in another study bacterial 16S rRNA 
genes amplified from two individual deep sea sponge samples were 
analysed by DGGE and the dominant bands were shown to be 
derived from thioautotrophic y-Proteobacteria [37]. These studies, 
together with the deeper pyrosequencing data presented here 
provide additional evidence for a common microbial community 
structure in deep sea sponges. 

AOA have previously been shown to be transcriptionally and 
metabolically active in nitrification in cold water sponges [16,29] 
and nitrification has also been detected in LMA sponges [42] , thus 
it seems likely that the AOA groups detected in these sponges may 
have similar metabolic capabilities. The other abundant group, 
thioautotrophic y-Proteobacteria has previously been found within 
sponges [11,37], and associated with other marine invertebrates 
such as marine tubeworms and sheMsh [38]. This class of bacteria 
appears to be a frequent symbiont, present in many marine 
invertebrates. Analysis of the seawater derived sequences revealed 
that at the 90% similarity threshold these y-Proteobacteria were 
essentially absent from seawater with a single read from one 
seawater sample grouping into the shared OTU. While it is possible 
that these bacteria are acquired from seawater they are clearly 
concentrated within the sponge tissue and appear to dominate the 
microbial communities of several sponge species. The SN sponge 
sample has a microbiota more typical of high-microbial abundance 
(HMA) sponges, with many bacterial species present and typical 
phyla such as Mtrospira, Acidobacteria, Poribacteria and Chloroflexi being 
observed. However this sponge also appears to have some 
similarities with the other deep water sponges under study; namely 
the presence of a common Tliaumarchaeota OTU and a Chromatiaks 
OTU, which together make up 19% of all reads in SN. 
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Figure 7. Relationship between microbial communities of deep 
sea sponges and Red Sea sponges by UniFrac UPGMA 
clustering. See Table 1 and Figure 6 for sample abbreviations for 
deep sea sponges. The pyrosequencing data for the Stelletta normani 
cortex (SNCo) and choanosome (SNCh) samples were analysed 
separately. Red Sea sponge samples SC-Sfy//ssa carteri, HE-Hyrtios 
erectus, XJ-Xestospongia testudinaria. Jackknife values are shown at 
nodes. The scale bar indicates the Unifrac distance. 
doi:1 0.1 371/journal.pone.0091 092.g007 



Evidence of a deep sea specific sponge community? 

The role of these Archaea and Bacteria within deep sea sponges is 
unknown, however the presence of two groups of microorganism 
whose closest studied relatives are chemolithoautotrophic is 
striking. The Thaumarchaeota are AOA, gaining energy through 
oxidation of ammonia while the sulphur oxidising y-Proteobacteria 
gain energy through oxidation of sulphide and other sulphur 
compounds. As previously mentioned, AOA have been shown to 
be major contributors to nitrification in cold water sponges and in 
the sponge Geodia harretti archaeal transcripts predicted to be 
involved in ammonia oxidation are highly abundant [16]. A 
symbiotic relationship between marine tubeworms and sulphur 
oxidising bacteria is well established, with tubeworms dependent 
on nutrients supplied by sulphur oxidising chemoautotrophic y- 
Proteobacteria [43]. A similar role for sulphur oxidising bacteria 
within sponges has not been shown although these bacteria appear 
to be widespread and abundant in deep sea sponges [37,38]. 
There is also a report of a specific putative symbiosis between free- 
living marine Thaumarchaeota and sulphur-oxidising y-Proteobacteria 
in sulphide rich mangrove swamps [44], implying that these two 
organisms themselves can form a close symbiotic relationship. 
How these common chemolithoautotrophic microorganisms 
contribute to the physiology of deep sea sponges is as yet unknown 
but their consistent presence impKes a potentially important 
symbiotic relationship. 



ai <" o 



5 2 .1] 

O g -D 



Conclusion 

This study presents an analysis of the shared prokaryotic 
diversity of four individual deep-water sponges and shows that the 
microbiota of these deep-water sponges share features with their 
shallow water counterparts. All sponge samples were found to 
contain diverse Bacteria and Archaea, and among the Archaea present 
a group classified as C. sjmbiosum was present in all samples. Three 
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of the individual sponges in the study L. diversichela, I. pellicula and 
P. compressa were classified as LMA sponges on the basis of their 
microbial communities. The microbial composition of the S. 
nonnani sample includes many OTUs from groups commonly 
found in HMA marine sponges. Pyrosequencing data indicated 
that the LMA sponges L. diversichela, I. pellicula and P. compressa were 
dominated by two microbial groups; AOA C. symbiosum group and 
sulphur-oxidising y-Proteobacteria, providing evidence of a putative 
common deep sea LMA sponge microbial community consisting 
of ammonia oiddising Archaea and sulphur oxidising y-Proteobacleria. 
The data presented here are derived from four individual sponge 
samples and additional data will be required to determine if this 
apparent deep sea LMA sponge microbiad community is shared 
with other deep sea sponges. Due to the nature of the current 
study the contribution of these common and apparently abundant 
members of the deep-sea sponge community to the shared 
metabolism of the sponge community is unknown but it is likely 
that these abundant microbes play an important role in sponge 
biology. 

Methods 

Sample collection 

Specific permission was not required, to obtain the marine 
sponge samples used in this study as they were collected in Irish 
territorial water, by an Irish research vessel, funded by the Irish 
government. The sponge samples do not involve endangered or 
protected sponge species. The sponge samples used in this study 
were collected with the remotely operated vehitie (R.O.V.) Holland 
/during the Biodiscovery cruise 2010 aboard R.V. Celtic Explorer. 
Upon retrieval the sponge samples were washed with sterile, 
artificial seawater (33.3 g/L Instant Ocean, Aquarium Systems ~ 
Blacksburg, VA, USA) and stored at — 80°G until molecular work 
was carried out in our laboratories in Cork. A part of each sample 
was also used for taxonomic identification by Christine Morrow, 
Queens University, Belfast. Additionally, a water sample was 
retrieved during a CTD (Conductivity-Temperature-Depth) mea- 
surement. 30 L water were collected as close as possible to the 
sponge samphng site, filtered through a 0.2 |xm membrane filter 
(Whatman - Austin, TX, USA) and the filter was immediately 
frozen at — 80°C. Depths and GPS location of samples are 
indicatc-d in Table 1 . Duplicate samples were not collected due to 
sampling difiiculties at such depths. 

Metagenomic DNA extraction from seawater and sponge 
samples 

DNA was extracted from filters using WaterMaster DNA 
Purification Kit (Epicentre Biotechnologies, Madison, WI, USA) 
according to the manufacturer's instructions. Extracted DNA was 
stored at -20°C. 

The sponge tissue (3-5 g) was cut into fine pieces with a sterile 
razorblade and then ground to a fine powder under hquid 
nitrogen using a mortar and pestie. For the S. normani sample, the 
cortex was first separated from the choanosome, cleaned carefully 
with sterile artificial seawater and any remaining sediment from 
the surface was removed with a sterile razor blade. For the 
choanosoin(; a cross section of the ball shaped tissue was taken in 
order to include inner- and outer areas of the choanosome. While 
processed separately reads from choanosome and cortex were 
ultimately pooled for some analyses. The ground sponge tissue was 
added to a lysis buffer adapted from Brady (100 mM Tris, 
100 mM EDTA, 1.5 M NaCl (w/v), 1% CTAB (w/v), 2% SDS 
(w/v); 5 ml buffer per 1 g sponge tissue; [45]) and incubated for 
2 h at 70°C. Metagenomic DNA was then extracted as previously 



described [1 1] . DNA solutions were analysed by gel electropho- 
resis, quantified by spectrophotometry (NanoDrop ND-1000- 
Wilmington, DE, USA) and then stored at -20°C. 

PCR amplihcon library preparation for pyrosequencing 

PCR amplicon hbraries of the V5-V6 region of 16S rRNA 
genes were prepared from all metagenomic DNAs. Universal 
primers U789f (5'-TAGATACCCSSGTAGTCC-3') and U1068r 
(5'-CTGACGRCRGCCATGC-3'), targeting both bacteria and 
archaea [13], were adapted for pyrosequencing by the addition of 
sequencing adapters and multiplex identifier (MID) sequences as 
per Table S 1 . Each 50 jtl PCR reaction comprised 1 x buffer, 
0.2 mM dNTPs (both Fermentas, St. Leon-Rot, Germany), 
0.1 |iM of each primer (Sigma Aldrich, Arklow, Ireland), 2 U 
Taq polymerase (Fermentas), ~10 ng template DNA and dH20. 
PCR cycle conditions were as reported previously [13]. To 
minimise PCR bias three individual reactions were performed per 
template and equimolar amounts of PCR products from each of 
the three reactions were pooled for pyrosequencing. PCR products 
were purified using Qiagen PCR Purification Kit (Qiagen Ltd., 
UK) as per the manufacturer's instructions. Barcoded samples 
were pooled and sequenced on a GS FLX Titanium platform (454 
Life S( i(;nc(;s) at the Uni\ ersity of Liverpool, Centre for Genomic 
Research, Liverpool, UK. 

Pyrosequencing data analysis 

Initial quality filtering and barcode assignment of reads was 
performed using the QIIME package [32] . Secjuences shorter than 
200 nucleotides, with an average quality score <25, or with >6 
ambiguous bases were removed from the analysis. AH reads were 
then analysed for errors and corrected using Acacial.52 [31]. 
Reads were then grouped into OTUs at 97%, 94% and 90% 
sequence similarities using uclust. Each OTU was then assigned to 
a taxonomic group using the RDP classifier at 50% confidence 
and the Greengenes database of assigned sequences [46] . OTUs 
were then aligned to the pre-ahgned Greengenes 16S data using 
PyNAST [47] . Potentially chimeric sequences were identified and 
removed from the dataset using ChimeraSlayer [48]. 

Rarefaction and rank abundance curves were calculated from 
OTU tabk-s gcrKTat(;d at 97'/o similarity using alpha diversity and 
rank abundance scripts within the QIIME pipeline. Shannon 
indices and Chaol species estimators were calculated from OTU 
tables generated at 90%, 94% and 97% similarity using alpha 
diversity scripts within the QIIME pipeline. 

Raw sequences are deposited in MG-RAST [49] and the 
NCBI short read archive with the following MG-RAST ID and 
BioSample numbers; IP - 4533583.3 SAMN02402455; LD 

- 4533578.3, SAMN02402456; PC - 4533577.3, SAMN02402457; 
SNCh - 4533579.3, SAMN02402458; SNCo - 4533580.3, 
SAMN02402459; W-1 - 4533582.3, SAMN02402460, W-2 

- 4533576.3, SAMN0240246 1 ; W-3 - 453358 1.3, SAMN02402462. 

Cloning and analyses of near full length 16S rRNA genes 

To generate longer sequences, suitable for more detailed 
phylogenetic analyses 16S rRNA genes were amplified from 
sponge metagenomic DNAs using primer sets 27f/1492r for 
eubacterial 16S rRNA genes [50] and Arch21F/Arch958R for 
archaeal 16S rRNA genes [51]. PCR amplicons were cloned into 
pJETl.2 (Thermo Scientific) and a selection of clones from each 
library were sequenced. Each of these libraries was found to 
contain sequences that were near-identical (>99%) to the most 
abundant reads in the V5-V6 libraries. Sequences were trimmed 
and aligned using PyNast [47] and analysed for chimaeras using 
ChimeraSlayer [48]. Phylogenetic trees were then calculated in 
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MEGA5 [52] using Neighbour-joining [53] and Maximum 
likelihood [54] algoridims. Near full lengdi I6S rRNA sequences 
of clones from this study have been deposited at GenBank with 
accession numbers Kr597097-KF597136. 

Supporting Information 

Figure SI Diversity of microbial communities in deep sea 
sponges and seawater. Rank abundance curve based on OTUs at 
97% similarity. See Table 1 for sample abbreviations. 

(TIF) 

Figure S2 Diversity of microbial communities in deep sea 
sponges and seawater. Rarefaction curve based on OTUs at 97% 
similarity. See Table 1 for sample abbreviations. 
(TIF) 
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